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Abstract. In the present paper the low density limit of the non-chronological 
multitime correlation functions of boson number type operators is investigated. We 
prove that the limiting truncated non-chronological correlation can be computed 
using only a sub-class of diagrams associated to non-crossing pair partitions and thus 
coincide with the non-truncated correlation functions of suitable free number operators. 
The independent in the limit subalgebras are found and the limiting statistics is 
investigated. In particular, it is found that the cumulants of certain elements coincide 
in the limit with the cumulants of the Poisson distribution. An explicit representation 
of the limiting correlation functions and thus of the limiting algebra is constructed in 
a special case through suitably defined quantum white noise operators. 
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1. INTRODUCTION 

The reduced dynamics of a quantum open system interacting with a reservoir in certain 
physical regimes is approximated by Markovian master equations. These regimes include 
the weak system-reservoir interactions and dilute reservoirs and in the theoretical 
framework they are described by certain limits. For a weakly interacting system one 
considers the limit as the coupling constant goes to zero (Weak Coupling Limit, WCL) 
whereas for a dilute reservoir one considers the limit as the density of the reservoir goes to 
zero (Low Density Limit, LDL) and an appropriate time rescaling should be performed 
in order to get a non-trivial limit. The Markovian reduced dynamics in these limits is 
considered in the review papers by Spohn and Lebowitz^l. The reduced dynamics in 
the LDL was considered in details later by Diimcke^ using the method based on the 
quantum Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy. 

The total dynamics in these limits is governed by various quantum stochastic 
equations. There is a unique up to now approach, called the stochastic limit method, 
which allows an efficient derivation of the stochastic equations in the WCL. This 
approach is based on the quantum white noise technique and was developed by Accardi, 
Lu, and Volovich 4 . 

The convergence of the evolution operator of the total system in the LDL to a 
solution of a quantum stochastic equation was proved by Accardi and Lu 5 and by 
Rudnicki, Alicki, and SadowskP. Recently the low density limit was investigated with 
the quantum white noise technique^'. This technique, well developed for the WCL, 
was non-trivially modified to include the LDL and for this case was called the stochastic 
golden rule for the low density limit. This technique was applied to the derivation of 
the quantum stochastic equations in the LDL. An advantage of the obtained equations 
is that they, in contrast with the exact Schrodinger equation, are explicitly solvable. At 
the same time they provide a good approximation of the exact dynamics. 

The approach of^uses the so called Fock-antiFock representation for the canonical 
commutation relations (CCR) algebra (this representation is unitary equivalent to the 
Gel'fand-Naimark-Segal representation). The difficulty with this approach is that the 
creation and annihilation operators in the Fock-antiFock Hilbert space do not describe 
creation and annihilation of physical particles and thus do not have direct physical 
meaning. To avoid this difficulty the investigation of the LDL directly in terms of 
the physical fields was performecP. Using this approach the chronological correlation 
functions in the LDL were found and the corresponding stochastic equations derived. 

In the present paper we investigate the low density limit of the non-chronologically 
ordered correlation functions of boson number type operators. The investigation is 
related with ab initio derivations of quantum stochastic equations describing quantum 
dynamics of a test particle interacting with a dilute gas. We find the limiting truncated 
correlation functions of the number type operators and show that they can be computed 
by representing the number operators through creation and annihilation operators and 
then considering only a sub-class of diagrams associated to non-crossing pair partitions. 
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This fact allows to represent the limiting truncated correlation functions as the non- 
truncated correlation functions of number operators of a free quantum white noise thus 
making a connection with the Voiculescu free probability theory. We find the limiting 
statistics and show that the cumulants of certain elements coincide in the limit with the 
cumulants of the Poisson distribution. 

The free probability theory was developed by Voiculescu around 1985 as a way 
to deal with von Neumann algebras of free groups. Then the theory was separated 
from this special context and began to develop as an independent field. In particular, 
applications of the free independence theory to random matrices were found. The details 
of free probability theory and its applications to random matrices could be found, for 
example, in reference d 10 ^ 11 ! 

Expectations of free random variables are characterized by diagrams associated 
to non-crossing pair partitions. The vanishing of crossing diagrams in the stochastic 
weak coupling limit for nonrelativistic QED and for the Anderson model was found 
in Refs^ ancP^, respectively, thus making a connection between the WCL and free 
probability. The WCL is typically described by the quantum Boltzmann statistics^. In 
Re£^ a generalized version of Boltzmann commutation relations, the so called entangled 
commutation relations, was found in the weak coupling limit for nonlinear interactions 
and possible applications to photon splitting cascades were discussed. 

The investigation of the multitime non-chronologically ordered correlation functions 
could have a connection with the behavior of fluctuations in certain asymptotic regimes. 
The latter is described in the review paper by Andries, Benatti, De Cock and Fanned. 
In that approach the limiting statistics is defined in terms of ground state distribution 
determined by non trivial pair partitions. The authors conjecture the appearance of 
exotic statistics in certain asymptotic regimes. The asymptotic fluctuations are the 
limiting correlation functions of appropriate centered elements and thus the results of 
the present paper could be applied to study the fluctuations in the low density limit. 

In Sec. II the truncated non-chronologically ordered correlation functions are 
defined and their low density limit is established (Theorem CD). In Sec. Ill the 
irreducible diagrams (pair partitions) which contribute to the limiting correlation 
functions are found (Theorem [2]). In Sec. IV the limiting truncated correlation functions 
are represented as correlation functions of a suitable free white noise. In Sec. V 
we identify the independent in the limit subalgebras (Theorem H]) and calculate the 
limiting cumulants which for some elements coincide with the cumulants of the Poisson 
distribution (Theorem[5]). In Sec. VI an explicit representation of the limiting correlation 
functions and thus of the limiting algebra is constructed for a special case by using 
suitable quantum white noise operators. 

2. THE CORRELATION FUNCTIONS IN THE LDL 

We begin this section with construction of a general class of non- commutative probability 
spaces relevant for the investigation of the low density limit. The framework of a *- 
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probability space is used. A relation between the objects defined in this section and the 
model of a test particle interacting with a dilute gas is given in Appendix A. 

Definition 1 A ^-probability space is a pair (A,uo), where A is a unital *-algebra 
over C and uj : A — > C is a state, i.e., a linear normalized, = 1, and strictly 

positive functional. 

Let H be a Hilbert space with inner product denoted by (•, •) (called as one 
particle Hilbert space), {S^gR a one parameter unitary group in Ti (a one particle free 
evolution), n a bounded positive operator in Ti (density operator) such that V£ G R, 
S- t hS t = h, and B a countable set of real numbers. 

Let r(7Y) be the symmetric Fock space over Ti. For any trace class self- adjoint 
operator T acting in Ti we denote by N(T) = dr(T) its second quantization operator 
in Y(Ti) and extend this definition by complex linearity to the set of all trace class 
operators T(Ti). For any T G T(Ti), uj G B, and a positive number e > we define the 
following operator in T(Ti): 

-itw/e 

AWO := ^—N(S t/£ TS_ t/£ ) (1) 

Let L(R) = H P eN L P (M), where L P (R) is the space of p-power intergable functions over R. 
For any open subset A C R let L(A) be the set of functions from L(R) with support in A. 
We denote by Aa, 6 the *-algebra generated by operators N TjU , j£ ((p) := J dt(p(t)N TjUlj£ (t) 
with T G T(H), uj G B, ip G 5(A) and denote A £ := Ar, £ . 

Let A^lg), g £ H be the creation and annihilation operators in T(7i) [we denote 
in the sequel A~(g) = A(g)} with the canonical commutation relations [A(f),A + (g)] = 
(f,g) and let ^4ccr be the algebra of polynomials in A ± (-). Any operator N(T) can be 
represented in terms of the creation and annihilation operators. For example, if T = 
|/) (si, where we use Dirac's notations for elements /, g G Ti, then N(T) = A + (f)A(g). 
An arbitrary operator N(T) can be expressed in terms of A ± using the fact that any 
trace class operator T is a limit of finite rank operators. Thus the algebra A £ is a 
subalgebra of ^4ccr- 

Let u)h be a gaussian gauge-invariant mean- zero state on ^4ccr with the two point 
correlation function Un(A + (f)A(g)) := (g,hf) (thus cJa(A r (T)) = Tr(nT) and here we 
use the assumption for T being trace class). Denoting by the same symbol its restriction 
to A\ t£ , we finally have for any e > and for any open subset A C R the *-probability 
space (^A, £ ,o; £ n)- 

Remark 1 The condition Wt: S^ t nS t = h leads to the invariance of the state t% under 
the free evolution generated by S t . 

With the notations above we define the non-chronologically ordered multitime 
correlation functions as 

W £ (ti,...,*„) ■■=u £ n{N Tl ^ £ (t 1 )...N Tn ^ £ (t n )) (2) 

W £ 

,h,Ti,uii,...,T n ,ui n 

(ipi, ...,(p n ):= u £ n{N TulJJU e{(Pi) ■ ■ ■ N TnjUntE ((p n )) (3) 
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We will use for the correlation functions (j2J) and ([3D also the shorter notations 
W e (ti, . . . , t n ) and W £ (<fi, . . . , (p n ). The reason for introducing the averaged operators 
Nt,u},e(¥>) and the averaged correlation functions ((3D is that, as we will show below, the 
non- averaged operators Nx tU j. £ {t) and the correlation functions ([2]) in the limit as e — > 
become singular distributions. Clearly, one has the relation 

W e (<p lt ...,cp n )= J dti . . . dt n W e (*i, . . . , . . . p n (tn) 

Definition 2 The truncated correlation functions W^(t\, . . . , t n ) are defined for n = 1 
fry H^(ti) := W £ (ti) and for n > 1 by induction through the relation: 

n 

w £ (t u ...,t n ) = wj{t u ..., tn )+j2 ZX T (^, • • • 



where Y^' is the sum over i\ < i 2 < . . . < ik 1} ifci+i < • • • < 2fc 2 , • • • ? *Jfcj+i < ■ ■ ■ <i n - 

The truncated correlation functions are often used in quantum field theory and in 
quantum kinetic theory^. They entirely determine the corresponding non-chronological 
correlation functions. Thus the investigation of the limit of the non-chronological 
correlation functions can be reduced to the investigation of the limit of the truncated 
correlation functions. 

We define the "projection" Pg := (27r) _1 J dtS t e~ ltE [it has the property PePe 1 = 
S(E — E')P E ] and for any k = 1, 2, ... ,n denote u k = u n + . . . + ojj,. The following 
theorem states the low density limit of the truncated correlation functions. 

Theorem 1 One has the limit in the sense of distributions in variables t±, . . . ,t n : 
Jim Wj^^^itu ■ ■ ■ , tn) = (2n) n - 1 5(t 2 -t x )... 5{t n - t n _ x ) 



J dE Tr [nP E+ a !1 T 1 P E+ a !2 T 2 . . . P £+tSn T n ] (4) 



x <5a 1( Q 

where Tr denotes trace and 8q u o is the Kronecker delta symbol. 
The theorem is a corollary of Theorem [2] from Section [31 

3. THE NON-TRIVIAL DIAGRAMS 

In the present section we investigate the low density limit of the non-chronologically 
ordered correlation functions for the particular case of operators of the form TJ = \fi)(gi\ 
and find the diagrams which are non-trivial in the low density limit. 

In order to simplify the notations we will use the following energy representation 
for the creation and annihilation operators: 

Aj := f^A + {P E M M ■= ^= £ A{S tl /e9i) 
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(a slightly different version of the energy representation was introduced irP). One has 
N Tu w u e{ti) = e~ itm / £ f dEiAfAi. Notice that the operator Aj is not the adjoint of A h 
The symbols Ai, A^~ are used only to simplify the notations below. 

A multitime correlation function can be expressed using Gaussianity of the state 
ujfi and the energy representation for the creation and annihilation operators as 

W e ,^Tx,wi,...,T n ,Wn(*l 5 ...,tn)= exp (-j ^ OOjti/ £ j ^ / dE\ . . .dEnLUeniA^AjJ . . . 

1=1 J 

x ^n(At k A hHn(A jk+1 At k+i ) . . . UJ £n (A jn A+) (5) 

where Yl' 1S the sum over k — 1, . . . , n, 1 = i\ < i 2 < ■ ■ ■ < ik, jk+i < ■ ■ ■ < jn, H < ji 
for I — 1, ... ,k and ji < i[ for I — k + 1, . . . , n. The sum contains terms of the form 

Ueh{K A h) ■ ■ ■ Uen(At k AlMn(A lk+1 At k+1 ) ■ ■ ■ ^en( A jn A t) (6) 

To each such term we associate a diagram by pairing in the string A± A1A2 A2 . . . A+A n 
the operators A^ and Aj t for I = 1, 2, . . . n. 

Definition 3 We say that the expression (0|) corresponds to a reducible diagram if there 
exists a nonempty subset I C {1, . . . ,n} (strict inclusion) such that %i E I ^ ji £ I. 
Otherwise we say that the expression (Ej) corresponds to an irreducible diagram. 

An important property of the truncated correlation functions (Def. [2]) is that they 
keep only all irreducible diagrams. The following are the examples of irreducible (first) 
and reducible (second) diagrams for n = 2: 

rm n n 

A+A 1 A+A 2 AtA 1 AiA 2 (7) 

Given an reducible diagram, one can represent the set {1, . . .n} as a union of several 
disjoint subsets I 1: . . . ,I t such that the diagram contains only pairings between operators 
with indices from the same subsets. In this sense a general reducible diagram can 
be represented as a union of mutually disjoint irreducible diagrams. Examples of the 
truncated correlation functions, the corresponding irreducible diagrams, and their limits 
as e — > for n — 1, 2, 3 are given below. 

Example 1 n = 1. The invariance of the state under the free evolution leads to the 
identity W?(t) = W £ (t) = W e (0) = ( 9l ,nfi). 

Example 2 n = 2. One has 

W^ihM) = W £ (t u t 2 ) - WefaWefo) = j dE 1 dE 2 uj £fi (AfA 2 )uj £h (A 1 A+) 

r e i(t2-ti)(# 2 -Ei)/e 

= j dE 1 dE 2 {g 2 ,P El hh)(g 1: {l + en)P E J 2 ) (8) 

This expression corresponds to the first (irreducible) diagram in ([7]) which is non-zero 
in the limit. Application of Lemma [T] (see Appendix B) to the r.h.s. of (jHJ) gives 

lim Wj(t u t 2 ) = 2vr5(t 2 - h) [ dE(g 2 , P E nfi) (gi, PeSz) 



cu efi (AtA 3 )uj eh (AiAt)uj E n(A2At 



Example 3 n = 3. One has 

W?{ti,t 2 ,t 3 )= JdExdE 2 dE 3 

+ uj £fl (A+A 2 )uj £h (A 1 A+)uj £ f l (A+A 3 ) 
This expression corresponds to the sum of the two irreducible diagrams: 



A+A 1 A+A 2 A+A 3 + A+A 1 A+A 2 A+A 3 
In this case only the first diagram is non-zero in the limit and Lemma Q] gives 

KmW?{t u t 2 ,t 3 ) = (2n) 2 5(t 3 -t 2 )6(t 2 -t 1 ) j dE(g 3) P E hh)( gi ,P E f 2 )(g 2 ,P E f 3 ) 
The case of arbitrary n is described by the following theorem. 

Theorem 2 Let Ti = \ fi){gi\, where fi, gi e H for I = 1, 2, . . . , n. One has the limit in 
the sense of distributions in variables t x , . . . , t n : 



lim W\ 



e->0 



xS, 



£,h,Tl,LtJ 



i,...,T„,w„(*l> ■■■>*») = ( 27r )™ <K*2 - *l) ■ ..5(t n - *„_i) 



2i,o / dE(g n , P E nfx)(gx, Pe+zJi) ■ ■ ■ (g 
For each n only the following irreducible diagram is non-zero as e — > 0: 



AtA 1 A+A 2 A+A 3 Ai . . . A n _ x A+ n A n 
Proof. Case (a): u x — ou 2 — . . . — ou n — 0. Using the correlation functions 

M^Afc.) = e^-^^ig^nfj 



(9) 



(10) 



^en(Aj A 



(g jg ,(l + en)f is ) 



one can write 



where 



as 



^ expji [(tx - t h )E x + ... + (t in - t jn )E tn ] /e] (e k F{E) + O 



(e k+l ))(ll) 



F(E) = l[{g jl ,P El nf il ) \{ (g jv Pe^) 

1=1 l=k+l 

Define the permutations pi and pj of the set (1, . . . , n) by pi(l) = i\ and pj(l) = ji 
for I = 1, . . . ,n and let p a = PiPj 1 - Consider the expression in the square brackets 
in the exponent in (II ip . The term proportional to t\ in this expression has the form 
ti(Ei — E ai ), where ai = p a (l). Thus ffTTl) can be written as 

^ expji [t n (E n -E an ) + ... + t^Ex - E ai j\ / e} (e k F(E) + 0(e k+1 ] 
and with the notations fli(E) = E n + . . . + E\ — E an — ... — E ai for I = 2, . . . , n as 



e k ~ 1 F(E) + 0(e k 



(12) 
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If the expression corresponds to an irreducible diagram then the functions Qi(E) 
are linearly independent and, since they are linear in their arguments, the convolution 
5(n 2 (E)) . . . S(Q n (E)) is well defined. 

In the case k > 1, since for any / = 2, . . . , n (see Lemma [1]): 

lim = 27r<y(t, - t,_i)£(n,(£)) (13) 

and — 1 > 0, the limit of (1T21) equals to zero. 

In the case k = 1 the expression ([6]) corresponds to the diagram (fTOl) and one has 

cj eft (,4+A l )u;^(A 1 A+) . . . cj efl (A n _ 1 A+) 



e i(t„-i„-i)n„(£)/e e i(t 2 -ii)n 2 (-B)/E 



(F(E) + 0(e)) (14) 



where Qi(E) — E\ — E\. Using (|T3|) one finds that the limit of the r.h.s. of (|T4|) is 

(27r)"- 1 5(t 2 -*!)■■■ <5(i n - ^05(^2 - £0 . . . 5(^ n - SO 

x PE 1 nfi)(g 1 , PeJi) ■ ■ ■ (g n -u PE n fn) 

Integration over E\. . .E n gives the equality ([9]) in the case (a). 

Case (b): arbitrary Ux, . . . , ui n . In this case the expression f fl~4|) in the 
decomposition ([5]) is multiplied by the factor exp(— i J2i uiU/e). The product can be 
written as 



e i(t n -t n -i){U n (E)-u> n )/e e i(t2-ti)(n 2 (E)-Cj2)/e 



e~ itiai/£ (F(E)+0(e] 



e e 
If a>i = then the statement of the theorem follows by the same arguments as in the 
case (a). If o>i 7^ then the limit of this term equals to zero by Riemann-Lebesgue 
lemma due to the presence of the rapidly oscillating factor exp(— iti^i/e). □ 



4. THE FREE WHITE NOISE NUMBER OPERATORS 

In the present section we show that the limiting truncated correlation functions coincide 
with the complete (i.e., non-truncated) correlation functions of the free white noise 
number operators. 

Definition 4 Free white noise operators Nf{t) are the operators satisfying the 
multiplication rule 

N T (t)N T ,(t') = S(t - t')N T , T r(t) (15) 

where the ^-product of any two operators T and T' is defined by T * T' : = 
2tt / dEP E TP E T'. 

Remark 2 We call the operators Nx(t) as free (or Boltzmann) number operators 

since they can be constructed using the creation and annihilation operators B^{t) 
satisfying the free relations Bj (t)B+(t') = 2n5(t — t')(f,g). In fact, define N\f^ g \(t) : = 
J dEBp E f(t)Bp Eg (t) and extend this definition by linearity to any T. Then such defined 
operators satisfy the relation (Th}) . 
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Let A be the algebra generated by the free white noise operators Nj-{t) and let 4>n 
be the state on A characterized by 0^(iVr(t)) = Tr (nT). 

Theorem 3 One has the equality 

lim^, Tl ,o,..,T„,o(^i> • • • ,tn) = <j>n{N T Ah) • • • N Tn (t n )) (16) 

Proof. By direct calculations using the Eq. (jl]) and the relation ffl5l) . 

The existence of the representation of the limiting truncated correlation functions 
by the free white noise number operators is related to the fact that only a sub-class of 
the non-crossing irreducible diagrams survives in the low density limit. We emphasize 
however, that the l.h.s. of Eq. ffl6|) is the limit of a truncated correlation function 
whereas the r.h.s. contains the complete correlation function. 

5. INDEPENDENCE AND THE GENERALIZED POISSON 
STATISTICS IN THE LDL 

The fact that the limiting truncated correlation functions are the distributions in 
variables tx, . . . , t n with support at ti — . . . — t n leads to the appearance of independent 
subalgebras in the low density limit. In the beginning of this section we remind the basic 
notions of independent subalgebras and of cumulants. Then we find the asymptotically 
independent subalgebras of A £ and discuss the limiting statistics. We show that the 
cumulants and the moments of certain elements in the algebra A e in the low density 
limit coincide with the cumulants and the moments of the Poisson distribution. 

Definition 5 Let (A, u) be a * -probability space. A family of unital ^--subalgebras 
{Ai}i(zi, Ai C A, is called independent if uj(ai . . . a n ) = whenever ai e A% 17 
ijj{a{) = 0, and k ^ I implies ik ^ k- 

Definition 6 Let (A, u) be a -^-probability space. Cumulants of the space (A, uS) are 
the multilinear functional n n : A n — ► C, n > 1, uniquely determined by Ki(a) : = 
u)(a), a G A, and for n > 1 by induction through the relation: 

u(a 1 ...a n )= 2j K \A j \({au---> a n)\A j ) 

7T,7T = :{Al,...,A fc } 

where the sum is over all partitions it of the set {l,...,n} and "(ax, . . . , a n )\A" 
designates the set of ai with i G A. 

Remark 3 The cumulants d$ for a ^-probability space (A £ ,u> £ n) are directly related to 
the truncated correlation functions. Namely, if a\ = N Tl)UJl>e {(pi), . . . , a n = N TntUj7it£ (ip n ) , 
then K { n\a u . . . , a n ) = W^fl, TliWlj ... iTniWn (yi, ...,(p n ). 

For the analysis of independence in the low density limit we introduce the notion 
of asymptotically independent subalgebras for a ^-probability space (A £ ,u £n ). 
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Definition 7 Let {A e ,oj e ^) be a * -probability space for the LDL. We say that a family 
of subalgebras A\ t£ , ■ ■ ■ ,Ai iE of A £ is asymptotically independent if 

\imu £ fy(ai, . . . ,a n ) = 

whenever ai G Ai u£ , u £fl (ai) = 0, and k ^ I implies ik ^ H- 

The next theorem identifies asymptotically independent subalgebras of A £ . 

Theorem 4 Let Ai, . . . , Aj be a family of disjoint open subsets in R. Then the family 
of subalgebras Aa- l , £ , ■ ■ ■ , A\ h£ is asymptotically independent. 

The proof follows from the fact that the truncated correlation functions become in the 
limit as e — > distributions in variables t%, . . . , t n with support at ti = t 2 = ■ ■ ■ = t n . □ 
Now let us analyze the statistics which appears in the low density limit. From 
Theorem [1] and the relation between the cumulants and the truncated correlation 
functions it follows that l-th cumulant for the element a = N^ui^if) i n the limit has 
the form 

lim^V,T,i^.--^) = ^«,o / dtdETrn[2n^t)P E T] 1 (17) 

e^O >!'''> 27T J 

We specify the further consideration to the case TC = L 2 (IR 3 ). Consider n = 1 
and St = e ltHl where Hi is the multiplication operator by the function u;(k) = 
|k| 2 , k G M 3 . Let T\ be an integral operator in 7i with the kernel T\(k, k') = 
(27r v /|k||k / |)~ 1 x [0 ^ ] (|k|)x [0iV ^ ] (|k / |), where A is a positive number and X[ ,^a] is the 
characteristic function of the interval [0, vA]. Let <fo(t) = (27r) _1 X[ 27r ](t)- 

Theorem 5 Let a\ = Nt x>u t£ {<Po)> where T\ and (f are defined as above. Then for any 
I G N one has 

Ki(a X} . . .,oa) = A(L,o 

or equivalently, the cumulants of the element a\ with uj = coincide in the low density 
limit with the cumulants of the Poisson distribution with expectation equal to A. 

Proof. The proof of the theorem is based on the direct calculation of the cumulants 
using Eq. (TT7T) . One has 
1 



Ki[a, ...,a) 



2tt 

One also has 



dt[27T<p Q (t)] 



j dETrn[P E T x } 1 = J dE J dk x . . . d^k^ 2 - E)^^, k 2 ) 

x 5(|k 2 | 2 - E)T x (k 2 , k 3 ) . . . <5(|k,| 2 - E)T x (k h kO 
= J dE[T x (VE,VE) J dk5(|k| 2 -E)]' = j dE XM (y/E) = A. 

Thus the r.h.s. of Eq. (TlTj) equals to one. This proves the theorem. □ 
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Moments of the element a\ with ui = in the low density limit are equal to the 
sum over all partitions of the limiting cumulants and given by Touchard polynomials: 

n 
k=l 

where S(n, k) is a Stirling number of the second kind, i.e., the number of partitions of 
a set of size n into k disjoint non-empty subsets. The limiting moments coincide with 
the moments of the Poisson distribution with expectation equal to A. For a\ one has 

\xrsxu) ef XaX) = B n 

where B n is the n-th Bell number, i.e., the number of partitions of a set of size n. The 
Bell numbers are the moments of the Poisson distribution with expectation equal to 1. 

6. AN OPERATOR REPRESENTATION OF THE LIMITING 
CORRELATION FUNCTIONS 

In the present section we explicitly realize the limiting correlation functions as 
correlation functions of certain operators acting in a suitable Hilbert space. Presence of 
delta functions in the limiting correlation functions suggests that they can be represented 
as correlation functions of certain white noise operators. Here such a representation is 
constructed in the special case using the results of^. 

Let go,g± G H satisfy the condition (g , S t g\) = for any t 6 R. Define for 
n,m = 0,1 the Hilbert space K nm := L 2 (Spec H\, d/i nm ), where SpecH 1 C R is the 
spectrum of H x and dfi nm := (g n , P E g n )(gm, PEng m )dE. Let K. := JC nm and let 

n, m=0,l 

1~~twN '■= r(L 2 (M, K)) be the symmetric Fock space over the Hilbert space of square 
integrable K- valued functions on R (abbreviation WN here stands for White Noise). 
Using the natural decomposition TCwn = <S> T(L 2 {R, JC nm )) one can define the 

n,m=0,l 

creation and annihilation operator valued distributions B^ n (E,t) acting in Hwn and 
satisfying the canonical commutation relations: 

[S- n (E, t), B+, n ,(E', t')} = 2n5(t' - t)8(E' - E) (g m , P E g m ,) (g n> , P E ng n ) (18) 

The operator valued distributions i?^ n (i?,t) are called time-energy quantum white 
noise due to the presence of 5(t' — t)5(E — E 1 ) in (fl8|) . Let define the number operators 

N m>n (E,t):= ln pf,a S B tn'( E ^) B n,n'( E ^) 

and denote N gm;9n (t) := J &E[N m , n {E,t) + B~jE,t) + B+ jTl (E,t)]. Let Q e H WN be 
the vacuum vector. 

Theorem 6 Let T x = \g rni ) (g m |, . . . , T k = \g mk )(g nk \, where m 1 ,n 1 , ... ,m k ,n k 6 
{0, 1}. One has the equality 

limW eATlA ...,T fc ,o(ti,...,4) = (V,N gmijgni (t 1 )...N gmk>gnk (t k )ty (19) 
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Proof, r.h.s. of (1191) has the form 

(^iV^^^i)...^^^) = j dE 1 ...dE k (Sl, [N muni (E 1 ,ti) + B~ umi (Ei,ti) 
+K uni (E 1 ,t 1 )] . . . [N mk>nk (E k ,t k ) + B- k!mk (E k ,t k ) + B+ k!nk (E k ,t k ))n) 

Let us denote N m , n (t) '■= J dEN mtn (E,t). The truncated correlation function 
corresponds to the term 

J dEdE' ( Q , 5~ mi (E,ti) N m2 N ma ,n 3 (£3 ) ■ • • N mk ^,n k Ah-i)B^ nk {E' , t k )Q) (20) 
Notice that N mjn (t)Q = 0. Therefore (1201 equals to 

J dEdE'(Sl, [. . . [[B~ imi (E,tx), N m2i n 9 (t 2 )], N m3ina (t 3 )] . ..N mk ^ nk ^t k ^)}B^ nk (E',t k )Q) 

The commutators can be calculated by induction using the canonical commutation 
relations (jTSj). The result is 

(2n) k - 2 5(t 2 -h)... - t fc _ 2 ) J dEdE'(g ni ,P E g m2 ) . . . (g 

x (n,B- k _ iimi (E,t k ^)B+ k)nh (E',t k )Q) (21) 

The last two-point correlation function can be calculated using the commutation 
relations (fISl) . This gives for (I2T1) the expression 

(2n) k - 1 5(t 2 -h)... 8{t k - t k -!) J dE(g ni ,P E g m2 ) ... (g )(g PEhg mi ) 

which coincides with the r.h.s. of ([9]) in the case w\ — . . . — to k — 0. □ 

Remark 4 The limiting correlation functions could be represented as expectations of 
certain quantum white noise operators in the general case if one could construct a 
Hilbert space Hwn, a vector £1 G Hwn, and operator valued distributions B^ g (E,t) 
and Nf t9 (E,t) in Hwn with the property Bj (E,t)ft = Nf t9 (E,t)Q = and satisfying 
the commutation relations 

[Bj g (E,t),B+ gl (E',t')} = 2n5(t' - t)6(E' - E)(f,P E f')(g',P E hg) (22) 

[Bj !g (E,t),N f , !g ,(E',t')} = 27rS(t'-t)S(E- E')(f,P E f)B glg (E,t) (23) 
[N Lg (E,t) : N fW (E' : t')] = 2n5(t' - t)6(E' - E)[(g,P E f)N ftg ,(E,t) 

-(g',P E f)N f ,, g (E,t)} (24) 

Suppose there exist such operators. Define Nf )3 (t) := / dE[Nf >g (E,t) + B~JE,t) + 
B~j (E,t)]. Then one can prove exactly in the same way as in Theorem^ that 

J™ W ^,ft,|/i)( ffl [,0,...,l/n>(9n[.0(*l, •••,*«) = • • • N fn,9n(tn)ty 
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APPENDIX A 

Here we make a connection between the objects defined in section II and the model of 
a test particle interacting with a dilute Bose gas (see RefPfor details). 

The one particle Hilbert space for this model has the form 7i = L 2 (IR 3 ), where M 3 
is the 3-dimensional coordinate or momentum space. The one particle free evolution 
is a unitary group St = e ltHl whose generator Hi in the momentum representation 
is the multiplication operator by the function u;(k) = |k| 2 /2m, where m is the mass 
of a gas particle. The test particle is characterized by its Hilbert space Hs and its 
free Hamiltonian H$ acting in Hs which is assumed to have a discrete spectrum. The 
discrete set B is the set of all transition frequencies of the test particle, or equivalently, 
the spectrum of its free Liouvillean — i[Hg, •]. 

The dynamics of a test particle interacting with a gas is described by an evolution 
operator U(t) acting in 7i s <g> T(H) and satisfying in the interaction picture, after the 
time rescaling t — > t/e, the following Schrodinger equation 



Here Qi tUJ are certain operators in Hs such that [Hs,Qi^] = —ujQi^ and Ti are certain 
operators in Ti. The explicit form of these operators is determined by the details of the 
microscopic interaction between the test particle and particles of the gas. Equation (1251) 
is the place where the operators N T ^^{t) appear. 

The condition S_ t nS t = n and positivity of n imply that for this model n is a 
multiplication operator by a function n : M 3 — > [0, oo). The value n(k) has the meaning 
of the density of gas particles at momentum k. If the state of the gas is u £fl then the 
density of gas particles and the rate of collisions between the test particle and the gas 
are of order e. Thus the limit e — > is the the low density limit. The limit is non-trivial 
since the dynamics is studied on the kinetic time scale of order \je. 




(25) 



APPENDIX B 



Let S(M.) be the Schwartz space over M and let S'(R) be the dual space of distributions. 
We reproduce the following lemma frorrP. 
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Lemma 1 One has the limit in S"(R) x S' 

pitx/e 

lim = 2Tr5(t)8(x) 

e~*0 e 

Proof. Let /, G S'(R) and let / be Fourier transform of /, f{r) = j dxe 1TX f(x). One 
has the identities 

/ e itx/e /■ /• /• 

dtdrr f(x)<f>(t) = lim / dr0(er) / dxe lrx /(a;) = lim / dr(f)(er)f(r) 
e e^Oj J 6^0 J 

Since / G S'(R), the function 0(er)/(r) satisfies the conditions of the Lebesgue lemma 
which allows to exchange the limit and integration in the last expression. Thus 

I = 0(0) f dr/(r) = 2tt0(O)/(O) □ 
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